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LIQUID  PROPELLANT  SMALL  ARMS 
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ABSTRACT 


This  is  the  final  report  on  Contract  No.  DAI-ll-022-0RD-P(8)  entitled 
"Interior  Ballistics  of  Liquid  Propellant  Small  Arms".  The  report  contains 
a  summary  of  a  literature  survey  of  the  liquid  propellant  field  and  the  re¬ 
sults  of  an  interior  ballistic  analysis  of  experimental  data. 

Experimental  data  obtained  by  Frankford  Arsenal  on  a  caliber  .60  weapon 
are  used  for  this  purpose  and  constants  for  use  in  the  ballistic  equations 
are  determined.  It  is  shown  that  a  quadratic  fora  function  and  a  linear  burn¬ 
ing  rate  equation  can  be  used  to  predict  the  maximum  pressure  and  muzzle 
velocity  with  sufficient  accuracy.  The  compromise  values  of  the  burning  rate 
and  form  factor  necessary  for  correlation  are  given.  Comparison  is  made  with 
solutions  obtained  with  a  cubic  form  function  equation  and  a  method  for  de¬ 
termining  these  form  factors  is  reported.  Data  from  rounds  with  various 
charge  to  projectile  mass  ratios  are  analyzed  end  show  that  the  burning  as 
indicated  by  a  fictitious  fractional  web  is  unaffected  by  a  charge  to  mass 
ratio  variation  of  2  to  1.  A  brief  analysis  is  given  indicating  the  apparent 
burning  surface  obtained  by  assuming  the  type  of  burning  rate  equation  used 
for  solid  propellants  but  using  a  value  of  n  >  1. 
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INTERIOR  BALLISTICS 


OF 

LIQUID  PROPELLANT  SMALL  ARMS 

I.  INTRODUCTION 

Research  on  the  use  of  liquid  propellants  and  the  development  of 
suitable  weapons  using  liquid  propellants  have  been  the  subject  of  pro¬ 
grams  sponsored  by  all  branches  of  the  service  for  the  past  several 
years.  These  programs  have  ranged  from  fundamental  studies  on  freezing 
points  of  various  LP  mixtures  to  the  development  of  40mm  repetitive  fire 
weapons.  Both  preloaded  and  injection  systems  have  been  investigated. 

The  research  activity  at  Frankford  Arsenal  has  been  directed  toward  the 
development  of  a  preloaded  monopropellant  system  using  hydrazine -hydrazine 
nitrate  in  a  caliber  .60  weapon.  Reproducible  performance  has  been  ob¬ 
tained  with  pyrotechnic  ignition  and  the  need  for  an  adequate  interior 
ballistic  system  has  been  expressed. 

The  Mechanics  Research  Department  of  the  American  Machine  &  Foundry 
Company  was  requested  by  the  Office,  Chief  of  Ordnance,  through  Frankford 
Arsenal,  to  conduct  a  research  program  leading  to  the  development  of  an 
interior  ballistic  theory  for  liquid  propellant  small  arms.  This  program 
was  begun  in  August  195^  on  Ordnance  Project  No.  TS1-11  under  Contract  No 
DAI-ll-022-0RD-(p)8 .  Work  continued  on  this  contract  until  its  termina¬ 
tion  November  30,  1955-  Two  simultaneous  phases  of  the  program  were  ini¬ 
tiated.  A  literature  survey  was  begun  to  keep  abreast  of  other  develop¬ 
ments  In  the  field  and  the  interior  ballistic  analysis  of  data  immediately 
available  from  Frankford  Arsenal  was  made .  The  information  obtained  from 

these  two  investigations  are  included  in  this  report. 
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II .  LITERATURE  SURVEY 


A  survey  of  existing  literature  served  to  acquaint  AMF  personnel 
with  the  status  of  LP  research  and  development  to-date,  and,  therefore, 
the  survey  was  continued  for  the  duration  of  the  contract.  The  survey 
was  to  serve  as  a  possible  source  of  ballistic  data  which  would  be 
analyzed  and  incorporated  in  the  interior  ballistic  theory  which  would 
be  developed  in  the  second  phase  of  the  project.  Unfortunately  much  of 
the  research  and  development  work  was  toward  the  development  of  operat¬ 
ing  guns  of  various  calibers  and  on  the  chemistry  of  the  various  liquid 
propellants.  Considerable  amounts  of  the  literature  concerned  caseless 
and  regenerative  feed  systems  which  are  of  only  cursory  interest  for  this 
project.  Thus  the  bibliography  in  Appendix  A  does  not  contain  all  of  the 
many  reports  which  were  surveyed,  but  only  those  which  contain  the  informa¬ 
tion  pertinent  to  the  interior  ballistics  of  pre-loaded  small  arms. 

A‘  Thermodynamic  and  Thermochemical  Properties  of  Hydrazine-Hvdrazl  ne 
Nitrate  ~  ” - — - * - 

1 •  Products  of  Decomposition 

It  is  generally  considered  that  hydrazine  completely  de¬ 
composes  into  and  Hg.  Any  water  present  acts  only  as  a 
diluent  and  does  not  participate  in  the  reaction.  The  reac¬ 
tion  when  hydrazine  nitrate  is  present  may  be  shown  as 
hN2H^  HNO^  ( aq)  — 5Ng  +  2N0  +  10HoC 

^4  +  2N0  - >-  2N2  +  2H20 

N2Hu - Ng  +  2Hg 
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These  reactions  assume  complete  dissociation  of  any  ammonia  which 
may  be  formed  as  an  intermediate  reaction.  Investigations  at 
Jet  Propulsion  Laboratory  (17)  have  assumed  that  ammonia  is  formed 
in  the  reaction  3^-^  4NH3  +  N^  3N2  +  6H.,,  and  have  shown 
that  the  maximum  impetus  is  obtained  vher.  approximately  35  per 
cent  ammonia  formation  occurs.  Per  cent  or  fraction  of  NH^  forma¬ 
tion  refers  to  the  fraction  of  the  theoretically  maximum  amount  of 
NH^  which  may  be  formed  by  the  above  reaction.  Thermodynamic  data 
are  presented  for  several  amounts  of  ammonia  formation.  Nozzle 
erosion  data  obtained  by  Armour  Research  Foundation  (3)  with  a 
vented  bomb  have  been  Used  to  indicate  gas  temperature  in  the 
throat  of  approximately  2000°K.  This  temperature  is  obtainable 
with  the  specific  composition  only  if  about  35  per  cent  ammonia 
formation  occurs.  These  data  give  direct  evidence  of  the  probable 
composition  of  the  gas  during  the  actual  ballistic  cycle. 

2.  Impetus  and  Isochoric  Flame  Temperature 

Values  for  the  impetus  and  isochoric  flame  temperature  for 
hydrazine  nitrate  water  have  been  reported  by  Amour  Research 
Foundation  (5>  6)  and  the  Jet  Propulsion  Laboratory  (9,,  17).  Armour 
Research  Foundation  assumed  no  dissociation  and  no  ammonia  formation. 

JPL  gave  data  with  up  to  2/3  of  the  theoretical  NH^  formation.  Im¬ 
petus  and  flame  temperature  increase  with  nitrate  addition  from  about 
31  x  10  ft-lb  per  lb  for  to  approximately  b6  x  10U  ft-lb  per  lb  at  85 

per  cent  NgHi-NO^.  Theoretical  flame  temperatures  vary  from  ll80°K 
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for  to  3375 °K  for  92.5$  N^H^NO^.  Both  impetus  and  flame 

temperatures  are  reduced  "by  the  addition  of  water.  It  has  been 
shown  (5,  17)  that  in  the  range  of  compositions  generally  used, 
l/3  ammonia  formation  increases  the  impetus  3  to  5  per  cent 
above  the  value  calculated  for  no  NH^  formation.  The  isochoric 
flame  temperature  is  more  sensitive,  increasing  approximately 
200 °K.  Hydrazine  alone  exhibits  a  maximum  impetus  with  about 
35  per  cent  NH^  formation  (17). 

3 •  Ratio  of  Specific  Heats  and  Co -volume 

The  ratio  of  specific  heats  is  also  dependent  upon  the 
products  of  combustion.  Values  for  Y  have  been  determined  by 
Armour  Research  Foundation  (6)  for  several  monopropellant 
compositions  ranging  from  100  per  cent  hydrazine  nitrate  to  100 
per  cent  hydrazine  and  for  several  different  fractions  of  water. 

In  these  calculations  it  was  assumed  that  no  ammonia  was  formed. 
The  values  of  7  vary  from  i.24l  for  O.Pmole  fraction  ~-2Eh°  HKO^ 
with  no  water  present  to  1.381  for  100  per  cent  N^H^.  The  pre¬ 
sence  of  ammonia  would  tend  to  lower  the  values. 

Several  values  of  Y  have  been  reported  in  literature.  Jet 
Propulsion  Laboratory  (l8)  used  Y  =  1.115  in  calculations  for 
a  3.2  in.  gun  using  a  mcnopropellant  of  71.7  per  cent  N2H4  and 
28.3  per  cent  N^H^NO^ •  This  value  seems  extremely  low  and  at 
best,  would  require  that  the  products  of  combustion  be  almost 
entirely  NH^  and  H^O  in  order  to  obtain  this  low  value.  Kerr- 
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Smith  (46)  uses  a  value  of  v  -  1.20  for  a  liquid  mixture  consist- 
lng  of  71  per  cent  hydrazine  nitrate,  13.6  per  cent  H2<),  6.7  per 
cent  N2Hv  and  8.5  per  cent  N02  and  products  of  only  H20,  and 
H2*  Armour  Research  Foundation  calculations  would  yield  I.34 
and  1.245  for  approximately  the  same  mixtures.  Adams  (l)  reported 
1.20  for  mixtures  containing  54  per  cent  hydrazine  and  37  -  4l  per 
cent  hydrazine  nitrate,  assuming  1  mole  NH^  formation  for  each 
00le  N2H4  not  oxidized  by  the  nitrate. 

The  co-volume  of  the  propellant  gases  is  also  a  strong  function 
of  their  composition.  Jet  Propulsion  Laboratory  (16,  18)  and  Armour 
Research  Foundation  (5,  6)  use  a  theoretical  value  of  38  in3  per  lb. 
for  low  pressure  equilibrium  with  large  quantities  of  NH^  formation. 
Kerr-Smith  (46)  uses  21.5  for  a  gas  mixture  assuming  no  ammonia  pre¬ 
sent.  Amour  Research  Foundation  (3)  also  reports  experimentally 
determined  values  of  30  for  a  35  per  cent  HNC>3  mixture  and  35 
for  a  95  per  cent  NgH^,  5  per  cent  H20  mixture.  These  values  were 
determined  in  a  closed  bomb  at  pressures  between  10,000  and  30,000 
psi.  At  this  pressure  level  the  effect  of  co-volume  is  small  and 
appreciable  error  in  the  determination  may  result.  Therefore,  these 
values  should  only  serve  as  a  guide. 

Linear  Burning  Rate 

Determination  of  the  linear  burning  or  consumption  rate  of 
hydrazine  nitrate  has  been  carried  out  by  several  investigators 
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including  Adams  and  stocks  (2)  in  England  and  the  United  states 
Naval  Ordnance  Laboratory  (to,  U)  and  Redel  Inc.  (St,  tt)  i„ 
this  country.  The  method  used  is  essentially  standard  and  con¬ 
sists  of  measuring  the  time  it  takes  the  liquid  in  a  glass  tube 
to  bum  a  given  distance  as  denoted  by  the  melting  or  breaking 
of  fine  vires  stretched  across  the  tube.  Tubes  of  3™  to  7™, 
have  been  used  but  the  most  extensive  vork  has  been  done  in 
5.0  -  5.5m  tubes.  An  increase  in  consumption  rate  with  tube 
diameter  has  been  attributed  to  variation  of  reaction  area  due 
to  the  change  in  meniscus  with  tube  diameter.  The  investigators 
have  found  that  there  is  a  critical  pressure  region  in  which  the 
propellant  will  not  ignite  or  sustain  turning .  The  level  end 
spread  of  region  depends  on  the  composition  of  the  propellant. 

The  upper  value  is  8500  psi  for  a  50  -  50  composition  hydrasine- 
hydrasine  nitrate  solution  and  increases  to  12,000  -  12,500  for 
a  25  per  cent  hydr.sine  nitrate  mixture.  In  general  the  consump¬ 
tion  rate  increases  with  nitrate  content  and  pressure.  For  each 

composition  then  has  been  found  at  least  one  region  of  negative 

slope  when  the  rate  sharnlv 

pxy  QnCreaseg  with  pressure.  Once  this 

region  is  past,  the  rate  again  increases  but  not  necessarily  at 
the  same  slope  as  before.  Typical  curves  reported  by  Redel  Inc. 
are  shown  in  Fig.  1.  The  slope  of  the  curve  in  the  log  log  pre¬ 
sentation  denotes  the  exponent  of  the  pressure  for  the  consumption 
rate  equations.  This  varies  from  approximately  2  to  nearly  3 
except  for  the  higher  nitrate  mixtures  at  pressures  of  15. 000  psi  or 
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FIG.  1  BURNING  RATE  AS  A  FUNCTION 
OF  PRESSURE  AND  COMPOSITION  (REDEL,  INC.) 
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higher.  A  slope  of  1  to  1.25  ls  shorn  for  these. 


Navord  Lab  (4l)  reports  consumption  rates  of  about  5  in  per 
ae=  for  pressures  around  4000  p.i  using  a  5-  tube  end  25  per  cent 

nitrate.  These  rates  are  extremely  high  compared  to  Redel- s  sort 

showing  5  ln/sec  at  25,000  to  nor, 

^  to  30,000  psl  for  approximately  the 

same  composition. 

Recently  limited  tests  (30)  have  shown  that  increased  sur¬ 
face  are.  caused  by  using  several  capillar  tubes  together  can 
double  the  consumption  rate. 

a  Effect  of  Stor«p-» 

Redel  Inc.  (29,  30)  reports  decreased  consumption  rates 
for  their  standard  composition  upon  storage.  A  decrease  of 
to  per  cent  in  3  months  was  reported  for  a  ample  containing 
aluminum  strip.  The  control  sample  decreased  25  per  cent 
but  the  sample  containing  stainless  steel  strip  showed  no 
decrease  until  after  six  month,  when  a  30  per  cent  decrease 

recorded;  at  this  time  the  control  ample  would  not 
ignite  at  30,000  psi. 

b-  Effect  of  Viscosity 

It  has  been  shown  (30.  to)  that  the  viscosity  effects  the 

consumption  rate.  Redel  Inc.  was  able  to  eliminate  one  refion 

of  negative  slope  from  an  established  burning  rate  curve  by 

increasing  the.  viscosity  with  the  addition  of  locust  bean  gum 
■fco  the  mixture . 
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c-  Effect  of  Initial  Temperature 

No  consumption  rate  determinations  have  been  reported 
at  reduced  temperatures.  Experimental  firing  data  at  re¬ 
duced  temperature  such  as  reported  by  Olin  Mathieson  (47) 
show  that  a  pressure  increase  occurs  below  about  30 °F  and 
reaches  a  maximum  at  the  freezing  point  of  the  mixture. 

Below  the  freezing  point  the  pressures  decrease,  but  be¬ 
come  erratic.  Pressures  increase  linearly  with  temperatures 
between  30°F  and  l40°F.  How  much  of  this  effect  is  due  to 
a  change  in  burning  rate  is  unknown. 

d.  Effect  of  Composition 

Increased  nitrate  content  increases  consumption  rate. 
Water  addition  tends  to  decrease  the  rate. 

5 •  Optimum  Composition 

Olin-Mathieson  Chemical  Corp.  conducted  a  firing  program  with 
a  20mm  weapon  to  determine  an  optimum  mixture  for  their  conditions 
It  is  reported  (12)  that  7.5  per  cent  water  and  20  -  25  per  cent 

N2H4-HN03  appear  to  give  the  maximum  velocity  to  pressure  ratio 
for  that  weapon. 

•  Performance  Characteristics 
1.  General 

The  general  shape  of  the  pressure-time  record  obtained  from  a 
high  loading  density  liquid  propellant  gun  differs  markedly  from 
the  curve  produced  by  the  burning  of  the  typical  solid  propellant 
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charge.  The  rate  of  pressure  rise  is  much  more  rapid  than 
that  of  an  equivalent  solid  propellant.  The  region  of  maxi¬ 
mum  pressure  then  takes  on  one  of  two  general  characteristic 
shapes.  Ideally  it  flattens  out  into  a  plateau  shaped  maxi¬ 
mum  before  decaying  monotonically  to  zero.  The  second 
characteristic  curve  is  a  "saddle"  shape  consisting  of  two 
peaks.  The  first  of  these  is  referred  to  as  the  "ignition 
peak"  and  comes  at  the  end  of  the  sharp  rise.  A  sharp  pres¬ 
sure  decrease  forms  a  "valley"  before  the  second  peak  after 
which  the  pressure  decreases  monotonically  as  in  the  ideal 
case.  The  cause  of  the  peaks  is  not  fully  understood  but  are 
controllable  to  some  extent.  In  general,  an  average  of  the 
two  peaks  results  in  a  pressure  close  to  what  might  be  obtained 
on  a  "plateau"  type  round.  They  can  both  be  of  about  the  same 
magnitude  or  if  the  initial  peak  is  high,  the  secondary  peak  is 
low  and  vice  versa,  ibis  phenomenon  is  quite  sensitive  to 

ignition  and  somewhat  affected  by  composition  and/or  consumption 

rate. 

Typical  of  pressure  records  obtained  from  LP  charges  even  in 
closed  bomb  studies  is  the  presence  cf  "hash"  or  high  frequency 
vibrations  on  the  record  which  may  be  of  the  sane  order  of 
magnitude  as  the  pressure  signal  itself.  Lew  loading  density  guns 
give  pressure  traces  similar  to  closed  bomb  traces.  This  might 
be  anticipated  since  it  is  essentially  a  constant  volume  burning 
with  negligible  projectile  travel  prior  to  all  burnt. 


AMERICAN  MACHINE  &  FOUNDRY  COMPANY 

- -  -  10  - 

CONFIDENTIAL 


Olin-Mathieson  Chemical  Corp.  (12)  m  order  to  obtsln  evlaence 
on  the  behavior  of  a  low  loading  density  charge,  performed  a  large 
number  of  firings  using  an  8cc  tests  chamber  with  different  sited 
barrels  attached,  and  with  different  weight  projectiles,  but  with 
the  charge  limited  to  a  maximum  of  2„0cc. 

For  hydrazine  monopropellant  with  5  per  cent  water  and  5  per 
cent  nitric  acid,  it  was  observed  that  for  each  charge,  regardless 
of  bore  diameter  or  projectile  weight,  that  the  maximum  pressure 
remained  constant.  This  effect  was  reported  consistent  with  the 
ballistic  process  attributed  to  the  low  loading  density  liquid 
monoprapeUant  charge;  .mmely  that  maximum  pressure  occurs  before 
appreciable  projectile  movement  takes  place.  Telocity  was  re¬ 
ported  to  vary  in  accordance  with  charge  and  projectile  mass. 
Poudrlm-tUh,  at  the  U.  S.  ..aval  Proving  Ground  has  reported  experi¬ 
mental  data  indicating  a  minimum  pressure  occurring  at  about  10  per 
cent  free  volume  Por  a  given  charge,  pressure  is  highest  at  high 
values  Of  free  volume  (lov  loading  density). 

3-  Effect  of  Additives 

Tabular  data  were  reported  on  the  effects  of  additives  by 
Olin-Mathieson  Chemical  Cory.  {12)  and  Syracuse  University  (W). 
Olin-Mathieson  Corp.  attempted  to  reduce  the  hash  on  the  pressure 
tracer  by  the  use  of  additives,  of  the  compounds  studied,  zinc 
nitrate  and  ferrous  chloride  additions  up  to  1  per  cent  resulted 
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in  hash  reduction  without  velocity  degradation.  Chamber  pres¬ 
sure  was  increased  appreciably,  however,  according  to  the  re¬ 
port. 

Syracuse  University  studied  the  effect  of  additives  on 
freezing  point  depression. 

4.  Effect  of  Viscosity 

In  order  to  test  the  effect  of  propellant  viscosity  on  gun 
performance,  Detroit  Controls  Corporation  (7)>  added  small 
quantities  of  locust  bean  gum  (up  to  1 . 5  per  cent  by  weight)  to 
the  propellant,  increasing  its  viscosity  from  2.9  to  a  maximum 
of  7200  poises.  While  there  was  no  improvement  In  gun  performance, 
the  change  in  viscosity  of  the  propellant  had  a  marked  effect  on 
the  combustion  process.  The  most  apparent  change  is  that  ignition 
has  been  adversely  effected,  and  erratic  combtistion  has  followed. 

In  many  of  the  firings,  particularly  those  firing  the  most  viscous 
propellants,  much  of  the  propellant  was  not  burned  and  remained  in 
the  barrel  and  combustion  chamber  after  firing. 

See  also  previous  section  on  linear  burning  rate. 

5-  Effect  of  Chamber  Geometry 

Both  Frankford  Arsenal  (42)  and  Reaction  Motors  (4o)  have  re¬ 
ported  work  on  chamber  geometry  effec'.  s.  The  Frankford  Arsenal 
work  is  discussed  in  Section  BbaBb  on  ignition.  Reaction  Motors 
report  results  from  different  chamber  configurations  and  special 
grids  for  promoting  turbulence  after  ignition  It  is  shown  that 
chamber  geometry  has  a  definite  effect  on  pressure  and  velocity. 
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60  — fect  of  Ignition  System 

a.  pyrotechnic  System 

The  experimental  studies  on  ignition  of  Haul* 

e  013  01  -Liquid  mcnopropel- 

lant  by  pyrotechnic  means  were 

ns  were  directed  at  optimizing  the  igni¬ 
tion  system  from  the  following  considerations: 

1-  Venting,  location  and  size  of  vents 
2.  Penetration 
3 •  Ignition  Pressure 

When  a  pyrotechnic  is  used  as  a  primer,  it  is  usuaUy 
located  in  a  tub.  with  either  radial  and/or  axial  splt  hole;. 

The  propellant  surrounds  the  primer  tube  and  ls  ignited  from  the 
flames  issuing  Prom  tb.  spit  hole  due  to  the  burning  of  the  pry„. 
technic.  The  ignition  of  the  priming  material  is  usually  ac¬ 
complished  by  uaing  a  primer  cap  containing  a  vely  small  amount 

Of  pyrotechnic  which  is  sensitized 

sensitized  by  either  a  mechanical  blow 

or  by  an  electrical  method. 

(2  )  Venting  Studies 

(4>  »gber  of  Radial 

Frankford  Arsenal  (33)  reported  that  the  use  of 
optimally  spaced  multiple  radlai  injection  regions  yield 
nore  plateau  shaped  p-t  curves  than  were  obtained  with 
single  regions.  In  addition,  it  makes  possible  the 

achievement  of  velocities  comparable  to  those  obtained 
with  axial  systems. 
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(b)  —  Ct  of  Vajyin g  Radial  Ignition  Orifice  Area 

It  was  observed  by  Frankford  Arsenal  (33)  that  the 
peak  pressure  and  velocities  decreased  as  the  area  was 
increased  with  the  optimum  p-t  shape  occurring  in  the 
group  having  the  greatest  orifice  area,  it  was  also, 
noted,  that  the  peak  pressures  in  all  cases  occurred  in 
the  ignition  position  of  the  curve. 

^  Effect  of  Radial  Vent  Locations 

Olin-Mathieson  Chemical  Corp.  (12)  studied  the  effect 
of  radial  igniter  flame  on  ignition  and  the  overall  bal- 
Ms-tic  characteristics.  A  standard  3cc  test  chamber  was 
used  with  a  150  grain,  caliber  -.30  projectile.  A  con¬ 
stant  length  primer  flash  tube  was  used  extending  the 
entire  length  of  the  chamber.  With  this  configuration, 
the  location  of  the  igniter  flame  could  be  varied  by 
varying  the  position  of  the  radial  flash  holes  along  the 
tube,  without  affecting  the  volume  or  configuration  of 
the  gun  chamber  or  the  tube  itself. 

In  all  tests  conducted,  the  maximum  chamber  pressure 
increased  as  the  position  of  the  igniter  flame  was  moved 
from  the  rear  to  the  forward  section  of  the  chamber. 

The  pressure-time  records  change  from  an  approximate 
plateau  shaped  curve  with  rear  end  ignition  to  a  sharply 
peaked  curve  with  front  end  ignition. 
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■d*  gffgct  of  pOSltion  of  Axial  Tuition  Orifice  in 
the  Fuel  Colnwr.  - - - - - - 

Frankford  Arsenal  (33)  reported  that  initially  it 
had  been  observed  that  the  ignition  peak  pressure  de¬ 
creased  as  the  region  of  axial  ignition  was  shifted 
forward  in  the  fuel  column.  The  results  reported  show 
that  contrary  to  expectations,  the  intermediate  length 
yielded  the  highest  ignition  pressure  and  the  optimum 
pressure-time  shape. 

In  cylindrical  combustion  chambers,  as  used  in 
this  work,  axial  ignition  systems,  in  which  the  primer 
gases  are  vented  through  a  hole  drilled  in  the  end  of 
the  flash  tube,  are  more  efficient  than  radial  systems 
in  which  primer  gases  are  vented  through  holes  drilled 
in  the  flash  tube  wall.  The  qualification  on  chamber 
geometry  is  inserted  since  it  is  considered  possible 
that  the  difference  is  due  tc  loss  of  ignition  energy 
to  the  chamber  walls  in  some  of  the  radial  tests 

I 

( 2)  Penetration 

(  - ,  gffgct  of  Primer  Vent  Area  to  Cavity  Volume  Ratio, 
v'  a 

In  the  case  of  Redelinc.  Z-l  pyrotechnic  reported 
by  Redelinc.,  (31)  penetration  increases  rapidly  with 
increasing  ratios  up  to  0.1.  Beyond  this  value,  it  does 
not  change  appreciably  except  for  a  pronounced  dip  at 
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0.224*  Tills  point  was  checked  by  a  repeat  firing  and 
found  to  be  valid.  With  FFG,  black  powder,  on  the 
other  band,  penetration  increases  less  rapidly  than 
for  Z-l  and  is  not  appreciable  effected  by  increasing 
values  of  the  ratio  above  0.4. 

Effect  of  Density  of  Loading,  /X  , 

Redel  Inc.,  (31)  reported  that  the  displacement  of 
the  primer  products  is  not  dependent  upon,  £  ,  when 
S  /V  is  held  constant,  although  penetration  was  least 
for  all  loading  densities  at  lowest  values  of  S  /v 

V 

It  can  be  svmmarized  from  (31)  .for  the  weights  of 
primer  material®  Redel  In<  .  investigated,  that  the 
maximum  Jet  penetration  during  0.3  to  0.5  msec  is  from 
1.4  to  1.75  inches  for  Redel  Z-l  and  0  9  to  1.35  inches 
for  FFG.  The  optimum  S,v/V  value  is  about  0,. 5  in"'1  and 
the  optimum  &  ,  is  about  0,5 

(3)  Ignition  Pressure 

Jet  Propulsion  Laboratory  (ll)  stated  that  one  of  the  main 
problems  in  the  ignition  of  a  monoliiuid  propellant  has  been  to 
limit  the  ignltior  p>.u#&ure  to  t.  t^nKCuUdu-s  high  value.,  In 
general,  the  ignition  pressure  hag.  'tended  to  be  tor  high . 

( a )  Effect  of  Additives 

Jet  propulsion  Laboratory  ( 1.1 )  found  that  for  &.  con¬ 
stant  charge,  varying  the  amount  >  f  RFNA  or  WFBA  igniter 
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is  not  a  satisfactory  means  of  limiting  the  ignition, 
pressure.  However,  by  adding  water  to  the  hydrazine 
the  ignition  pressure  could  be  limited  in  a  narrow 
range  of  water  concentration. 

(h )  Effect  of  Chamber  Geometry 

Frankf Ox  d  Arsenal  {^2^  observed  while  conducting 
tests  with  cased  arid  caseless  rounds,  that  the  ignition 
peaks  were  consistently  much  higher  for  the  caselees 
than  comparable  eased  firings.  It  was  postulated  that 
two  basic  differences  might  contribute  to  increased 
ignition  pressure  encountered.  The  first  of  these 
lies  in  the  slightly  larger  diameter  of  the  easeless 
chambers  which  might  permit  more  propellant  to  be 
ignited  by  the  initial  primer  exhaust  per  unit  length 
of  penetration.  A  second  lies  in  the  difference  in 
the  radii  of  curvature  of  the  rear  of  the  ease  and 
of  the  caseless  breech  plug  obturate.  The  curvatures 
are  such  that  the  amount  of  propellant  to  the  roar  of 
the  axial  orifice  in  the  caseless  breech  plug  is  2.1cc, 
whereas  the  corresponding  volume  in  curved  configuration 
is  6.1cc. 

(4)  Pressure  Variations 

The  U.  S.  Naval  Proving  Ground  at  Dahlgren  (lj>)  postulated, 

from  experimental  observations,  that  the  quantity  of  energy 
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spiled  to  the  propellant  by  the  Igniter,  the  mf  of 

deUTe,T  °f  thiS  «»  ^~on  of  thX.  energy 

to  the  propellant,  and  the  Meant  of  turbulence  created 

in  the  propellant  by  the  Igniter  must  he  controlled  vitVr 

narrow  limits  If  satisfactory  unifo,«ty  Is  to  he  obtained. 

^5)  Perchlorate  Ignition 

Olin-Mathieson  Chemical  Corp.  (lit)  folmi  that  thev 
could  obtain  the  beat  Ignition  of  38  grans  of  liquid  propel.. 
innt,  by  ualng  1/2  gran  ammonium  perchlorate  W!en  perchlorate 
la  need  as  a  booster  In  the  gin.  It  Is  partisan,  deposed 

^  thi“  Prl”r  fl“8  “  **  resulting  nurture  of  orygen  rich 
flene  pin.  decomposed  perchlorate  is  blasted  into  the  hydras 
in  the  vicinity  of  the  priaer. 

Theoretical  Considering 

Bedel  Inc,,  (32)  derived  at  the  following  relation  to 
describe  the  dynamics  of  a  pyrotechnic  jet; 

v.  g(piV  “* 

K 


where : 


.  /(Pl-Pe,)  K 

V - —  t 

'  2fL  Cc 


Ks  *  coefficient  of  total  *  / 

steady  flow  at  the  Inst^Lt^aVSc^y 

V  *  velocity  of  jet 

C  "  unsteady  flow  constant 
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L  = 

P  - 

"t  ss 


fluid  density- 
conduit  test  length 
pressure 
time 


b-  Spark  Ignition  System 

Considerable  effort  has  been  expended  in  the  field,  not 
only  m  developing  a  suitable  spark  ignition  system,  but 
in  theoretical  analyses  of  the  mechanism  of  sparking  and 
ignition  from  sparx  sources . 

Theoretical  Considerations 

Stanford  Research  Institute  (*)  is  conducting  a 
liquid  propellant  program  in  which  they  investigating 
the  mechanism  by  which  a  spark  ignites  the  monopropellant . 

The  framework  of  their  hypothesis  is  provided  by  two 
experimental  observations,  that  the  proximity  of  the  walls 
to  the  spark  makes  ignition  easier,  and  that  the  shock 
waves  are  provided  by  a  very  high  energy  spark.  The  rate 
of  chemical  reaction  caused  by  the  spark  will  have  an  in¬ 
fluence  on  the  behavior  of  the  outward  moving  shock  wave 
by  virtue  of  its  energy  contribution  to  the  expanding  piston. 
The  reaction  tone  may  travel  with  the  head  of  the  shock  wave, 
since  the  shock  strength  is  so  great  that  the  temperature  in 
the  shock  front  is  appreciable.  However,  unless  a  steady 
state  detonation  wave  is  set  up,  the  spherical  shock  wav. 
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will  attenuate  as  It  moves  outward  and  the  reaction  will 
proceed  at  a  slower  rate  behind  the  shock  front  in  the 
region  of  elevated  pressure.  The  behavior  of  the  shock 
wave  will  be  determined  by  the  rate  of  input  of  energy 
to  the  piston  by  the  spark  and  by  the  chemical  reaction. 

It  is  supposed  that  if  the  pressure  in  the  neighborhood 
of  the  reaction  falls  below  the  minimum  pressure  for 
burning,  the  reaction  will  cease.  It  is  further  supposed 
that  if  rigid  walls  are  situated  at  a  suitable  distance 
from  the  spark,  the  increased  pressure  in  the  first  re¬ 
flected  shock  will  sustain  the  pressure  in  the  neighbor¬ 
hood  of  the  reaction  zone  for  a  time  sufficiently  long  that 
the  critical  volume  of  liquid  can  have  an  opportunity  to 
react. 

The  hydrodynamic  problem  then  is  one  of  the  calcula¬ 
tion  of  the  behavior  of  a  spherical  shock  wave  generated 
by  a  spherical  piston  moving  outward  at  a  rate  determined 
by  the  rate  of  delivery  of  spark  energy  and  the  rate  of 
reaction  of  the  propellant. 

Redel  Inc.,  (27)  hypothesized  that  the  mechanism  of 
electrical  action  consists  of  the  following  three  steps: 

1.  heating  and  vaporization  of  a  portion  of  the 
liquid, 

2.  ignition  of  the  vapor  in  the  arc, 

3.  additional  energy  used  to  promote  turbulence 
and  ionization. 
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As  soon  as  the  switch  is  closed,  a  potential  is 
established  across  the  spark  gap.  As  the  solution  is 
heated,  the  resistance  of  the  liquid  decreases ,  and  the 
current  flow  Increases. 


At  some  time  the  path  between  the  electrodes  con¬ 
tains  a  quantity  of  vapor  such  that  the  electrical  flow 
changes  in  character  from  a  flow  of  electrons  through  a 
liquid  to  a  flow  of  ions  through  a  gas  phase,  constitut¬ 
ing  an  arc.  This  can  occur  only  if  the  voltage  at  this 
instant  is  still  above  minimum  spark  voltage  for  that 
gap  and  gas  constant.  As  this  arcing  process  is  essen¬ 
tially  a  high  temperature  one,  a  portion  of  the  propellant 
vapor  is  raised  from  the  boiling  temperature  to  the  decomposi¬ 
tion  temperature.  An  unknown  quantity  of  the  available 
energy  is  employed  in  this  step  and  is  believed  to  be  very 
small. 


After  the  arc  is  formed,  the  energy  remaining  in  the 
system  is  poured  through  the  arc.  This  continues  until 
the  voltages  at  the  electrodes  drops  below  the  value  re¬ 
quired  to  maintain  an  arc  across  the  gap.  The  arc  then 
dies  leaving  a  small  residual  of  energy  in  the  system  The 
characteristics  of  the  electrical  system  determine  the 


amount  of  residual  energy.  The  energy  being  put  into  the 
arc  is  used  in  the  third  step  of  the  mechanism  in  heating 
of  the  reaction  product  gases,  dissociating  product  gases 


at  high  temperature,  bringing  additional  propellant  vapor 
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to  decomposition  temperature  and  even  vaporizing  additional 
propellant.  These  action*  are  all  exhibited  as  the  crea¬ 
tion  of  turbulence  around  the  spark  gap. 

Redel  further  states  (20)  that  experimental  data 
point  to  the  conclusion  that  the  greater  part  of  the  spark 
discharge  energy  is  in  the  formative  phase  and/or  line 
loss.  It  is  felt,  however,  that  the  majority  of  energy 
in  the  spark  phase  is  available  to  promote  turbulence  in 
the  system.  Also  a  comparison  of  pyrotechnic  and  spark 
Ignition,  as  applied  to  monopropellant  initiation  in¬ 
dicates  that  the  type  of  Ignition  Is  not  important  so  long 
as  a  certain  minimum  energy  is  reached. 

University  of  Michigan  (37)  found  that  when  an  impulse 
potential  is  applied  to  a  set  of  electrodes  submerged  in  a 
conducting  fluid,  the  current  discharge  takes  place  in  two 
distinct  phases.  The  first  phase,  called  the  "fornative 
phase”,  is  characterized  by  a  high  voltage  drop  across  the 
electrodes  and  a  low  current  flow.  The  second  phase  called 
the  "arc  phase",  la  characterized  by  a  low  voltage  drop 
across  the  electrodes  and  a  high  current. 

The  energy  delivered  to  the  monoliquid  propellant  dur¬ 
ing  the  arc  phase  of  discharge  does  not  seem  to  be  the  sole 
criterion  for  ignition.  Pour  ignitions  were  obtained  with 
arc  energies  between  4.5  and  8  Joules,  and  yet  no  ignition 
was  obtained  on  the  other  tests  when  the  arc  energy  was  8  to 

AMERICAN  MACHINE  A  FOUNDRY  COMPANY 

CONFIDENTIAL  -  22  - 


f 


CONFIDENTIAL 


10  Joules  energy.  Thus  it  appears  that  the  formative 
energy  and  total  spark  energy  dissipated,  as  well  as  the 
arc  energy,  all  bear  some  relation  to  the  ignition  process. 

(2)  Factors  Affecting  Spark  Ignition 

It  was  found  by  Olin-Matfaieson  Chemical  Corp.  (10,  12), 
that  ignition  of  hydrazine  would  not  occur  when  a  spark 
was  discharged  beneath  the  liquid  surface,  even  with  the 
chamber  completely  filled  with  liquid  propellant  and 
sealed  by  a  heavy  aluminum  diaphram.  Ignition  occurred 
only  when  the  chamber  was  prepressurized  with  nitrogen  to 
approximately  2000  psi .  (See  section  on  electrodes). 

(a)  Plenums 

Hydrazine  and  hydrazine -hydrazine  nitrate  mono¬ 
propellant  mixtures  have  been  ignited  successfully 
by  means  of  a  single  spark  discharge  without  prepres¬ 
surization  by  simply  confining  a  small  quantity  of 
propellant  in  a  plenum  chamber  surrounding  the 
electrode.  A  spark  plug  was  used  incorporating  a 
plenum  volume,  and  the  0.30  caliber  test  gun  was  fired 
with  the  system  using  the  propellant  mixture  containing 
60  per  cent  hydrazine,  33  per  cent  hydrazine  nitrate 
and  7  per  cent  water.  The  minimum  spark  energy  required 
for  ignition  was  6  Joules  (3  mfd  condenser  to  2000 
volts) . 
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(b)  Formative  Energy 

According  to  the  University  of  Michigan  (36),  it 
seems  fairly  definite  that  successful  ignition  can  be 
obtained  with  an  all  formative  discharge.  Twelve  gun 
firing  tests  were  made  with  all  formative  and  10  of 
these  resulted  in  successful  ignition.  The  minimum 
energy  for  ignition  seems  to  be  slightly  below  15 
joules  for  the  particular  system  being  used,  which 
is  the  same  minimum  energy  as  for  the  same  standard 
formative -arc  discharge.  The  peak  power  delivered 
during  the  all  formative  is  in  the  same  range  as  that 
obtained  in  the  formative-arc  discharge  since  the 
latter  has  its  peak  during  the  formative  period. 

(c)  Electrodes 

Univeristy  of  Michigan  (38)  showed,  in  a  10  ml 
bomb,  completely  filled  with  hydrazine -hydrazine -nitrate 
liquid  propellant,  that  in  order  to  obtain  ignition  it 
is  necessary  to  cover  all  but  the  tips  of  the  electrodes 
with  an  insulating  material  to  prevent  an  undue  amount 
of  stored  energy  from  being  dissipated  by  conduction  dur¬ 
ing  the  foriative  stage. 

A  particular  difference  between  blunt  and  pointed 
electrodes  was  the  amount  of  scatter  in  the  energy  data 
(35)-  For  any  given  spark  gap,  the  energy  obtained  with 
blunt  electrodes  was  considerably  more  erratic,  parti¬ 


cularly  in  the  formative  phase. 
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In  comparing  the  changes  in  formative  energy  with 
other  properties  of  the  propellant,  it  was  found  that 
electrical  conductivity,  viscosity,  and  formative 
energy  are  all  affected  proportionally  by  changes  in 
temperature  (39).  By  utilizing  various  strengths  of 
potassium  chloride  solutions  and  maintaining  a  con¬ 
stant  propellant  viscosity,  it  has  been  shown  that  an 
exponentional  relation  exists  between  formative  energy 
and  conductivity.  This  same  sort  of  relationship  holds 

true  for  hydrazine  propellant,  but  the  exponent  is  not 
the  same. 


It  has  been  observed,  in  spark  discharge  tests 
made  with  blunt  and  pointed  electrodes  that  both  energy 
and  duration  appear  to  increase  slightly  with  increas¬ 
ing  viscosity  (39). 

(e)  Voltage 

Initially  it  was  determined  (8)  that  the  release  of 
approximately  100  joules  of  electrical  energy  from  a 
condenser  bank  was  sufficient  for  ignition.  Tests  were 
made  at  both  2000  and  4000  volts  to  determine  the 
effects  of  capacitor  voltage  on  round  performance  while 
maintaining  constant  energy  in  the  capacitors.  No 
appreciable  change  in  performance  was  observed  other  than 
a  slight  increase  in  peak  combustion  chamber  pressure, 
accompanied  by  a  minor  improvement  in  reproducibility,  in 
tests  made  at  4000  volts. 
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C.  Interior  Ballistics 


It  is  necessary,  in  any  interior  ballistic  analysis,  to  consider  the 
following  five  relations: 

1 -  Energy  Balance 
2.  Equation  ol'  State 
3-  Equation  of  Motion 
^  •  Propellant  Burning  l,»w 
5 »  Form  Function 

The  first  three  relations  are  essentially  the  same  for  any  type  analysis. 
However ,  items  4  and  5,  the  burning  law  and  for.  function,  dan  dlffer  con. 

siderably . 

JPL  (16)  presented  a  method  of  Interior  ballistics  for  use  on  a  high 
speed  calculator  (REfiC)  In  order  to  correlate  information  on 

1-  e?fert1ofUSo^2“'aCteriS';Ira  including 

burning  raS°''  PreSS,"e  <“1  e"ecti,n! 

2.  form  function,  the  mode  of  propellant  break  up  and 
resulting  change  of  burning  rate, 

3*  effective  impetus, 

effective  loss  to  barrel 

5-  pressure  dependence  to  burning  rate. 

Heat  loss  was  considered  to  be  proportional  ir  ,  ..  ,  , 

~  -  vjjui  Lioiua,. l  to  une  kinetic  energy 

Of  the  projectile,  as  accounted  for  in  solid  propellant  interior  ballistics. 

The  burning  law  used,  was  of  the  type  commonly  applied  to  the  burning 

of  solid  propellants,  the  fraction  nf  woi, 

'  ctuon  on  web  unourned  is  proportional  to 

pressure  raised  to  a  power,  as 


D  ~/V' 

at 
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where; 


o(  =  pressure  exponent 

P  *  "breech  pressure  (assumes  combustion  takes  place  at  breech 

pressure) 

J3  =  burning  rate  coefficient 

D  =  web  size 

f  =  fraction  of  web  unburned 

Since  liquid  propellant  has  no  web,  in  the  solid  propellant  sense, 
the  quantities  f,/2  and  D  were  redefined  in  terms  of  the  physical  pro¬ 
perties  of  the  liquid.  Thus,  D  was  defined  as  the  quotient  of  the 
liquid  voxume  divided  by  the  surface  area  exposed  to  the  igniter.  Then 
if  is  the  usual  burning  rate  coefficient,  D(df/dt)  will  be  the  linear 
burning  rate  normal  to  the  surface,  at  least  initially.  Since  D  can  be 
obtained  empirically  through  the  ballistics,/*  can  therefore  be  determined. 

It  was  postulated,  that  in  order  to  obtain  a  flat  pressure  curve  in  a 
liquid  gun  it  is  necessary  to  have  initially  regressive  burning  followed 
by  progressive.  The  fora  function  used  to  obtain  this  condition  was 
taken  as 

(p  -  (1-f)  (l+Of)  (1+Ojf) 

where : 

(p  ■  fraction  of  propellant  burned 

0  and  0^  are  constants  and  00^>0  for  the  conditions  mentioned. 

After  the  propellant,  has  eomplot  °ly  burned,  the  expansion  was  considered 
as  adiabatic. 
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Computations  were  made  using  a  quadratic  form  function  (0^  a  0) 
and  assuming  that  burning  occurs  at  space  mean  pressure.  These  com¬ 
putations  showed  that  dependence  of  the  various  ballistic  parameters 
with  variations  in  each  other  and  other  dimensionless  parameters. 
However,  the  results  were  too  preliminary  to  draw  any  significant  con¬ 
clusions  . 

Redel  (28)  presents  a  method  whereby  the  experimental  pressure 
time  curve  is  analyzed  to  determine  the  relation  between  propellant 
burned,  propellant  burning  surface  area,  temperature,  velocity  and 
travel  and  time.  Although  this  method  may  present  Insight  into  the 
values  of  the  various  ballistic  parameters,  at  the  conditions  at  which 
the  experimental  data  was  taken,  in  its  present  form  it  cannot  be  used 
as  a  method  for  weapon  design. 

The  method  presented,  does  utilize  the  necessary  relations  of 
energy  balance,  equation  of  motion  and  equation  of  state.  However, 
unconventional  interior  ballistic  notation  is  used. 

The  total  heat  loss  was  evaluated  from  the  muzzle  conditions,  and 
then  were  inserted  linearly  with  time  The  internal  energy  was  assumed 
to  be  either  a  linear  or  a  quadratic  function  of  temperature.  An 
equation  of  state  is  used  in  which  gas  imperfections  are  accounted  for 
by  a  compressibility  factor,  Z ,  rather  than  the  usual  covolume  correc¬ 
tion. 
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JPL  (lS)  also  reports  a  method  for  calculating  a  system  assuming  a 
linear  pressure  rise,  a  constant,  pressure  period  at  the  maximum  until 
°1-1  c:c*  :~L-  expansion.  Equations  are  presented  for 

each  of  four  periods  in  the  ballistic  cycle.  This  has  been  achieved 
by  forcing  the  desired  pressure  function  into  the  ballistic  equations, 
and  actual  pressure  functions  only  approach  the  desired  conditions  at 
best. 

Kedel  Inc.,  (29)  reports  a  more  radical  treatment  of  the  burning 
by  studying  the  propellant  break-up.  Use  is  made  of  Taylor's  stability 
theory  to  determine  the  minimum  and  most  probable  droplet  sizes.  An 
attempt  is  made  to  take  into  account  the  surface  tension  and  viscosity 
effects.  Reacting  surface  areas  of  60  times  the  chamber  cross-section 
are  calculated. 

HI-  DEVELOPMENT  OF  INTERIOR  BALLISTIC  THEORY 
A.  Basic  Equations 

An  interior  ballistic  calculation  for  either  a  solid  or  a  liquid 
propellant  weapon  involves  the  simultaneous  solution  of  four  basic 
equations .  ' 

The  energy  balance  equation 


NRT 

o 


7~  1 


NRT 

r-  1 


? 


(1) 
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The  Abel  equation  of  state 


-*]- 


NET 


(2) 


The  equation  of  motion 


AP  =  m." 


The  burning  rate  equation 


dt  "  f s(B/2)pn 


dv 

dt 


(3) 


W 


The  first  three  equations,  well  known  In  solid  propellant  halllstlos, 
cm  be  applied  directly  to  liquid  propellants  with  little  or  no  nonfic¬ 
tion.  The  burning  rate  equation  is  based  on  the  interface  burning 
surface,  s  ,  vhlch  will  vary  during  burning.  With  solid  propel  lets, 
there  is  a  constrained  geometrle  surface  and  a  definite  relation 
exists  between  the  surface  exposed  to  burning  and  the  fraction  of  the 
propellant  .a,  consumed.  This  relationship  is  generally  referred  to  a, 
the  foe  function.  Kith  liquid  propellants,  there  is  no  surface  constraint 
and  the  interfaclal  surface  my  vary  considerably  without  regard  to  the 
fraction  of  propellant  mass  consumed.  Thus  the  evaluation  of  a  suitable 

expression  for  the  burning  equation  consltutes  the  maior  problem  in  the 
liquid  propellant  ballistic  theory. 

For  solid  propellants  the  form  function  is  expressed  as  a  function  of 
th.  fractional  wab,  f,  which  is  the  ratio  of  the  charaoterlstlc  weh  length 
of  the  grain  to  its  initial  value,  for  liquid  propellants  the  fractions! 

•weh  can  have  no  such  physical  significance.  However,  it  is  convenient  to 


AMERICAN  MACHINE  & 


FOUNDRY  COMPANY 


CONFIDENTIAL 


30 


consider  it  ss  a  continuous  variable  which  Is  Initially  "1»  and 

mono  tonic  ally  decreases  tc  zero  at  "«.-n  vr  +  >. 

°  at  a"1  burnt  •  Using  this  notation 


and  the  Eq.  (4)  becomes 


S  "  “v 


(5) 


-W  —  -  Rpn 

W  dt  ~  BP 


If  It  1,  determined  that  the  form  function  varies  with  ballistic 
parameters  evaluation  of  s  must  be  performed  from  Eq.  (k) . 

The  progresslvity  of  the  burning  is  determined  by  the  change  of 
the  interfacial  surface  and  can  be  easily  determined  by  (5). 

The  literature  survey  has  shown  that  the  following  factors  may 
affect  the  burning  rate; 

1.  chamber  geometry 

2.  loading  density 
3 •  ignition  system 

propellant  properties 
5-  expansion  rate 

Reaction  Motors  Inc, ,  (4q)  ha*  * 

'  '  9>  *  Showa  that  Pre“ignition  agitation  of 

the  propellant  does  not  affect  the  bailie  tics . 

B-  Working  Equations 

In  developing  ballistic  equation,  for  solid  propellants,  it  has 
been  the  practice  to  asstmre  that  the  pressure  exponent,  n,  in  the  burn¬ 
ing  equation  is  unity.  This  simplification  is  justified  for  solid  pro- 
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pellants  since  n  varies  between  G„7  and  1.0  for  most  propellants.  Data 
on  liquid  propellants  have  indicated  that  n  may  be  between  2  and  3  for 
the  hydrazine-hydrazine  nitrate  monopropellant  compositions  generally 
used.  However,  it  was  believed,  that  in  the  interests  of  simplicity 
an  interior  ballistic  theory  based  on  a  linear  burning  rate  equation 
should  be  attempted.  Therefore,  taking  n  =  1,  the  four  basic  equations 
are  combined  into  the  general  ballistic  equation  in  terms  of  the  fol¬ 
lowing  dimensionless  parameters ; 


i  BT 

I  *  V  CF  V*  dimensionless  velocity 
-  jl/k.  ,  dimensionless  travel 

"o*  dimensionless  pressure 


— j —  (B/W),  *  \l  CF/id"  1/Vb,  dimensionless  burning 

rate  constant 


t-  d/ac0)  ^ 


fcf. 
i  et  c 


dimensionless  time 


The  resulting  ballistic  equation  becomes 


where : 


(~k^)  I2  -  0  -  1/2  / 


~  -  aA  01 
P  o  a  A 


0  =  0  (f) 


f  =  f  -CM 
o 
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uj  _  <i{l2) 

T  d  >. 


d?  = 


a 


(y-D  = 


dl/^r  =  d  ^  /2I 

7  -  ^ 

d+^95  a  - 1) 

heat  loss  proportionality  factor. 


(8c) 

(8d) 

(8e) 

(8f) 


The  system  of  Eq.  (8)  may  be  solved  if  the  form  function  Eq.  (8a) 
is  given.  Solution  of  this  system  for  a  quadratic  form  function  in  f 
is  well  known.  It  will  be  shown  later  that  a  cubic  form  function  beBt 
describes  the  burning  for  the  plateau  shaped  pressure  time  curves,  al¬ 
though  experimental  data  can  be  closely  approximated  with  a  quadrat i c 
function  by  modifying  certain  other  parameters.  Solutions  for  both 
cases  are  reported. 

C •  Idealized  Pressure  Time  Curve 

Experimental  data  have  shown  two  fundamental  type s  of  pressure¬ 
time  relationship:  the  saddle  shape  and  the  plateau  shape.  It  has 
heen  assumed  that  the  plateau  type  curve  which  consists  of  a  nearly 
linear  pressure  rise  to  maximum  or  close  to  maximum  followed  by  a 
period  of  nearly  constant  pressure  Before  s  monotonic  decrease  during 
the  latter  part  of  the  cycle  is  the  desirable  uniform  burning  pattern, 
and  the  saddle  shape  in  which  two  peaks  are  obtained  constitutes  varia¬ 
tion  in  the  burning  due  to  the  effect  of  various  parameters.  The 
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plateau  type  of  curve  is  idealized  and  expressions  for  the  surface 
function  are  derived  in  Appendix  B,  assuming  a  linear  pressure  rise 
and  then  a  constant  pressure  until  all  burnt.  This  analysis  shows 
that  the  variation  of  surface  with  fractional  web,  f,  can  be  expressed 
by  two  straight  lines,  as  shown  in  Fig.  2,  representing  constant  pro- 
gressivity  in  each  region.  The  plateau  region  has  the  lower  pr ogres - 
sivity  as  would  be  expected.  The  two  solid  straight  lines  may  be 
closely  approximated  by  a  parabola  as  indicated  by  the  dotted  line. 
Based  on  conventional  formulation  of  form  functions,  any  stright  line 
in  Fig.  2  represents  a  quadratic  form  function  while  a  parabola  with 
the  directrix  oriented  parallel  to  the  f  -  axis  represents  a  cubic 
form  function  of  the  form 


0  =  k0  -  kif  +  V2  +  k3f3  (9) 

where:  kQ,  k^,  Is.^  and  are  constants. 

An  analysis  of  the  magnitude  of  these  form  factors  is  given  in 
Appendix  C.  It  is  shown  that 

=  1  (10a) 

=  (1  +  kg  +  k3)  >  0=  kx  at  f  =  0  (10b) 

-  1  -  k2  -  2k3  >  0  -  k5  at  f  -  1  (10c) 

Other  criteria  for  the  validity  of  the  form  factors  are  given  also. 
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FIG.  3  REPRESENTATIVE  PRESSURE -TIME  DATA 
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FIG.  5  EXPERIMENTAL  AND  CALCULATED 


PRESSURE -TRAVEL  CURVES 


For  the  calculations  and  in  analyzing  the  experimental  data 
the  Impetus  of  the  propellant  was  taken  as  334,000  ft  and  Tj  vas 
assumed  to  be  30  in^/lb. 

The  constant  kjin  the  quadratic  form  function  is  an  indica¬ 
tion  of  the  progressivity  of  the  propellant  since 

Sc*^f  U3) 

and 


The  modified  specific  heat  ratio ,T  ,  reflects  the  energy  in 
the  gases  since 


E  .  NT 

0  -n 


Therefore  a  high  value  of  /  means  that  there  is  less  energy  in 
the  generated  gas.  These  factors  are  apparent  in  Figs.  4,  5  and  6. 
An  increase  in  kg  tends  to  flatten  out  the  pressure -travel  area  and 
give  a  lover  peak  for  the  ssae  linear  burning  rate  .  The  increase 
in  y  causes  the  tail  end  of  the  curve  to  become  lover  and  is  parti¬ 
cularly  effective  after  all  burnt. 


A  comparison  of  the  (fi  -  V  relationship  obtained  from  experimental 
data  with  values  calculated  from  the  quadratic  form  function  is  shown 
in  Fig.  7*  The  individual  points  on  the  dashed  ^e  experimental 
values  while  the  solid  curves  represent  matching  data  from  a  quadratic 
form  function.  These  data  show  that  0  function  follows  a  quadratic  form 
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to  *  velocity  of  2000  fpa  where  the  experimental  data  shows  a 
pronounced  departure  from  the  quadratic  form  with  a  decrease  in 
slope.  The  matching  values  show  that  for  7  -  1.4,  kg  =  0.948 
and  0-=  0.804.  This  agrees  with  the  pressure-travel  curve  except 
for  the  value  of  kg. 

Analysis  of  the  experimental  data  using  ~Y  =  1-4  showed  that 
all  burnt  occurs  at  shot  ejection  and  0.65-  However,  it  is 
shown  that  good  agreement  between  experimental  and  theoretical 
curves  can  be  obtained  by  compromising  on  the  quadratic  form 
factor  and  assuming  a  higher  burning  rate  factor. 

2.  Cubic  Form  Function 

Ballistic  calculations  were  made  using  the  cubic  form  function 
as  expressed  by  Eq.  (9).  Values  for  kg  and  k^  were  determined 
from  experimental  data  using  the  criteria  set  forth  in  Eq.  (10) 
group.  Since  kr  represents  the  slope  of  the  surface  curve  at 
initial  conditions,  a.  value  of  0.05  was  assigned  for  the  first 
calculations.  The  results  of  these  calculations  are  shown  in 
Fig.  8.  These  calculations  show  that  the  shape  of  the  curve,  maxi¬ 
mum  pressure  and  area  under  the  curve  can  be  closely  approximated 
with 

k„  =  0.05 

kg  =  -0.528 

k3  -  0.739 
!T~  =  0.76 
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FIG.  8  COMPARISON  OF  CALUCLATED  AND  EXPERIMENTAL  PRESSURE- 
TRAVEL  DATA  USING  CUBIC  FORM  FUNCTION 
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The  error  in  maximum  pressure  is  less  than  4  per  cent  and  the 
muzzle  velocity  from  integration,  is  only  1.2  per  cent  higher. 

This  curve  is  shifted  in  time  from  the  experimental  curve. 

A  higher  value  for  the  initial  slope,  k^,  was  chosen  in  order  to 
increase  the  rate  of  pressure  rise.  The  initial  slope  for  the 
quadratic  form  was  determined  to  he  0.11  and  a  value  of  0.1  was 
chosen  for  the  next  set  of  calculations .  The  results  shown  in 
Fig.  9  indicate  that  this  results  in  a  sharper  peak  hut  that  the 
deviation  in  peak  pressure  is  small  and  a  compromise  on  values  of 
kg,  k^,  and  CT’can  keep  the  deviations  to  less  than  5  per  cent. 

This  data  also  uses  ^ =  0.69  which  is  close  to  the  experimentally 
determined  value  of  0-65 

Either  the  quadratic  form  function  or  the  cubic  equation  may 
he  used  for  calculating  the  pressure,  velocity  and  travel  for  a 
desired  weapon  based  on  the  performance  of  the  caliber  .60  weapon. 
Unfortunately,  data  available  on  other  weapons  were  insufficient  for 
good  analysis  and  comparison  purposes. 

A  comparison  of  the  7  -f  relation  as  determined  from  experimental 
data  with  data  calculated  by  quadratic  and  cubic  form  functions  is 
shown  in  Fig.  10.  The  cubic  form  function  shows  somewhat  better 
correlation  but  is  mere  difficult  to  use  in  the  ballistic  equations. 

If  electronic  computers  are  available,  the  added  complexity  of 
the  cubic  form  function  is  unimportant,  however,  it  does  involve  the 
evaluation  of  one  additional  constant  for  the  form  function.  Since 

for  engineering,  applications,  adequate  agreement  of  maximum  pressure, 
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FIG.  10  COMPARISON  OF  0  VS.  f 
FROM  EXPERIMENTAL  AND  CALCULATED  DATA 
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muzzle  velocity  and  travel  can  be  obtained  with  a  quadratic  func 
tion,  use  of  the  cubic  probably  cannot  be  justified  on  the  basis 
of  slightly  better  correlation. 


^ '  Effect  of  Charge  to  Projectile  Mass  Ratio 

Fraiikford  Arsenal  personnel  fired  a  series  of  rounds  in  which 
the  charge  was  held  constant  at  the  established  value  and  the 
projectile  mass  was  reduced,  yielding  a  variation  of  the  ratio 
of  charge  to  projectile  mass  of  as  much  as  2  to  1.  Representative 
pressure-time  data  for  each  C/m  are  shown  are  shown  in  Fig.  11. 
These  data  show  that  the  magnitude  of  the  ratio  has  negligible 
effect  on  the  maximum  pressure  or  the  general  shape  of  the  curve. 
The  smaller  projectile  mass  results  in  a  higher  velocity  and  hence 
a  shorter  time  to  shot  ejection  and  greater  expansion  rate.  These 
data  were  analyzed  also  on  the  basis  of  the  relationship  between 
0,  the  fraction  of  gas  formed  and  f,  the  fictitious  fraction  of 
the  web  consumed.  It  is  assumed  that  f  -  1  -  (v/vj.  These  data 
are  show  in  Fig.  12.  It  is  apparent  from  these  data  the  charge 
mass  ratic  and  ..he  resulting  changes  in  expansion  rate  (or  pro¬ 
jectile  velocity:  does  not  change  the  0  -  f  relationship. 

Discussion  of  Data  Evaluation 


Two  parameters  in  the  ballistic 


equations  are  not  well  defined 


by  experimental  data  to  date, 
ratio  of  specific  heats,  . 


These  are  the  co-volume,  r?  nnri  the 
Covolumes  reported  in  the  literature 


™ry  ““drably  from  a  mpxima,  of  36  ln3  per  lb  a  ^  of 
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in  per  lb.  depending  on  the  assumed  exposition  of  the  product. 

Of  combustion.  Using  38  i„3  pen  lb.  nested  In  extrxely  high 

calculated  pressure,  end  a  value  of  30  in3/lb  was  flMlly  8elected 

as  giving  reasonable  results.  This  value 

8  value  has  been  reported  from 

experiments  also. 

Fit.  literature  it  woSd  be  expected  that  7  should  be  between 
1.2  and  1.3.  An  attempt  to  determine  7  from  experiments  data 
fails  because  it  was  not  possible  to  obtain  sufficient  informa¬ 
tion  after  all  burnt.  In  the  ballistic  escalations  7-  ,  lA  mB 

used  because  it  gave  good  correlation  with  the  experiments  p-x 
data. 

The  first  set  of  data  analyzed  eave  . 

ya  gave  excellent  results.  Travel 

“  velocity  determined  by  integration  agreed  with  indicated  shot 

ejection  time  end  muscle  velocity.  The  second  set  of  data,  for 

the  rounds  fired  with  varying  charge  to  mass  ratios,  showed  discrep- 

•ncies.  Rounds  fired  with  the  sxe  charge  and  projectile  mas,  as 

the  previous  data  showed  longer  snot  ejection  times.  When  integrated, 

however,  the  data  indicated  shot  ejection  prior  to  the  signs  on  the  ’ 

pressure  trace  and  was  in  agreement  with  the  earlier  data.  It  was 

— ed,  therefore,  that  the  muscle  blanking  unit  was  not  recording 

properly  and  the  time  scale  obtained  from  integration  was  used  for 
all  rounds. 

The  "malfunction"  of  the  muscle  blank  made  it  possible  to  obtain 
limited  pressure  data  after  shot  ejection.  RvSuation  of  the  0 
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function  from  these  data  indicated  that  probably  is  close  to 
the  1.2  value  expected  as  indicated  in  Fig.  13.  However,  since 
a  value  of  1.4  must  be  used  in  the  theoretical  calculations  to 
match  experimental  pres sure -travel  data,  this  has  little  effect 
on  the  choice  of  form  function  constants  to  be  used  if  the 
assumption  of  delayed  muzzle  blanking  is  valid.  It  has  been 
impossible  to  prove  or  disprove  this  assumption.  Future  ex¬ 
perimental  work  should  be  carried  out  with  sufficient  tube 
length  to  insure  the  all  burnt  condition  prior  to  shot  ejection. 

analysis  OF  BURNING  SURFACE  FROM  EXPERIMENTAL  DATA 

If  Eq.  4  is  solved  for  S,  the  variation  of  S  with  time  which  would  give 
the  measured  pressure -time  curve  may  be  determined  assuming  this  burning  law. 

Remembering  that  N  =  C0,  Eq.  4  may  be  written 


d  0 


dt 


SB_  n 
C  r 


from  which 


R  _  SC  1  d0 
f  B  pn  dt 

Experimental  values  of  0  (t)  may  be  determined  from  Eq.  (ll). 


(15) 


(16) 


The  variation  of  0  thus  calculated  when  "7=  1.2  is  shown  in  Fig.  13. 

The  variation  of  S  which  would  be  required  to  give  the  measured  pressure  time 
curve  for  these  conditions  is  shown  in  Fig.  l4,  for  n  =  1  and  n  =  2.25.  The 
effect  of  assuming  n  =  2.25  is  to  flatten  the  first  peak  and  accentuate  the 
later  maximum.  ‘This  indicates  that  there  is  no  large  break-up  of  the  liquid 
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FIG.  13  VARIATION  OF  FRACTION  BURNED  WITH  TIME 
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FIG.  14  VARIATION  OF  RELATIVE  SURFACE  AREA  WITH  TIME 


during  the  early  burning  period.  The  large  increase  in  surface  area  occurs 
quite  suddenly  after  the  velocity  has  reached  approximately  2800  fps.  At 
this  point,  only  about  35  per  cent  of  the  charge  remains.  The  combination 
of  increased  chamber  volume,  smaller  liquid  charge,  projectile  velocity 
and  physical  properties  of  the  liquid  at  this  temperature  and  pressure  have 
apparently  reached  a  critical  point  allowing  break-up. 

Figure  1?  shows  a  typical  saddle-shaped  pressure-time  curve.  The  gas 
production  function,  <p  ,  and  the  surface  area  variation  are  shown  in  Figs. 

16  end  17.  The  surface  area  curve  differs  only  slightly  from  Fig.  13  which 
was  obtained  from  a  relatively  smooth  pressure  curve.  Since  both  surface 
curves  are  obtained  by  measuring  the  slope  of  an  experimentally  determined 
curve,  some  discrepancy  may  he  expected  from  the  reading  accuracy.  The  im¬ 
portant  information  revealed  by  these  data  is  the  trend  of  the  surface  area 
at  the  time  the  saddle  occurs  in  the  pressure.  In  the  time  0.6  to  1.0  msecs 
while  the  pressure  is  decreasing,  the  surface  area  is  actually  increasing 
appreciably.  The  first  surface  peak  corresponds  with  the  second  pressure 
peak,  however.  Thus  it  would  appear  that  the  saddle  shape  results  from  a 
change  in  the  linear  burning  rate  or  in  the  pressure  exponent  and  not  from 
a  sharp  decrease  in  burning  surface.  An  analysis  of  the  break-up  forces  would 
require  knowledge  of  the  viscosity,  surface  tension  and  density  of  the  liquid 
at  the  prevailing  temperature  and  pressure.  The  temperature  of  the  liquid 

would  be  a  function  of  the  thermal  properties  which  are  unknown:  estimates 
would  be  required. 
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v,  CONCLUSIONS 

The  development  of  an  interior  ballistic  theory  for  liquid  monopropellants 
is  by  no  means  conplete.  A  considerable  number  of  closely  controlled  experi¬ 
mental  firings  are  needed  to  provide  sufficient  data  on  the  effect  of  various 

parameters.  However,  several  conclusions  may  be  drawn  from  the  work  reported 
here. 

These  are: 

A.  The  assumption  of  a  linear  burning  rate  equation  with  either  a  qua¬ 
dratic  or  a  cubic  form  function  may  be  used  to  predict  the  pressure- 
travel  relationship.  For  the  quadratic  form  function,  k2  -  0.89  and 
0=  0.80.  Close  agreement  is  also  shown  for  the  cubic  form  function 
if  k5  -  0.05,  k2  -  -0.528,  k3  -  0.739  and  cr-  0.76.  When  k^  -  0.10 
good  agreement  can  be  obtained  with  k2  -  -0.508,  k3  -  o.70l*  and 

<7-»  0.69.  A  modified  specific  heat  ratio,  >-1.1*  must  be  used  in 
each  case. 

B.  The  quadratic  form  function  with  the  above  constants  and  co-volume 
assumption  provides  adequate  correlation  for  engineering  purposes. 

Even  with  the  benefit  of  electronic  computers,  the  use  of  a  cubic 
form  function  does  not  appear  justified  on  the  basis  of  these  data. 

C.  The  form  function  as  determined  from  the  experimental  data  at  high 
loading  densities  is  unaffected  by  the  ratio  of  charge  weight  to 
projectile  mass,  or  therefore,  the  expansion  ratio. 

D.  Analysis  of  the  burning  surface  by  Eq(l*)  indicates  that  under  the 
usual  burning  conditions  large  increases  in  the  burning  area  do  not 
occur  until  after  peak  pressure  has  been  reached. 
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S.  this  analysis  also  indicates  that  tha  saddle  shaped  variation  of  the 
experimental  pressore-time  data  may  be  due  to  a  sudden  decrease  in 
either  the  constant  or  the  pressure  exponent  in  the  burning  equation 
and  not  due  to  a  collapse  of  the  available  burning  surface. 

F.  More  accurate  determination  of  the  co-volume  is  necessaxy  to  ixprove 
the  ballistic  calculations.  The  pressure  is  very  sensitive  to  this 
parameter  for  weapons  using  high  loading  densities. 

Further  study  of  conclusions  D  and  E  should  be  made  using  ex¬ 
perimental  data  in  which  other  variations  in  the  pressure  time  data 
are  obtained.  If  the  analysis  is  valid,  it  should  be  determined  if 
the  sharp  increase  in  surface  area  after  peak  pressure  can  be  correlated 
with  some  hydrodynamic  and  flow  theories  which  have  be.  advanced. 

RECOMMENDATIONS 

Available  experimental  data  is  not  sufficiently  complete  „„  any  0ne 

system  to  result  in  an  interior  ballistic  theory.  Data  obtained  on  on.  system 

can  seldom  be  correlated  in  ballistic  equation,  with  data  from  another  system 

because  variations  in  several  parameters  generally  occur,  m  order  to  fully 

evaluate  the  effects  of  chamber  configuration,  initial  temperature  viscosity, 

ignition,  loading  density,  and  charge  to  projectile  mass  ratio  in  the  ballistic 

equations,  experimental  data  must  be  obtained  from  a  series  of  closely  controlled 

experiments.  The  experimental  weapon  must  be  assigned  so  that  each  of  the 

parameters  may  be  varied  without  changing  any  of  the  other  parameters  mentioned 
above . 

The  following  specific  recommendations  are  made? 

A.  An  experiments!  evaluation  of  the  co-volume  for  the  range  of  mixtures 
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generally  used.  The  evaluation  should  be  made  at  pressures  In 
excess  of  50,000  psi  in  a  closed  bomb. 

B.  Analytical  and  experimental  evaluation  of  the  heat  losses  during  the 
ballistic  cycle  should  be  made. 

C.  Conduct  a  closely  controlled  experimental  program  in  which  each  of 
the  factors  suspected  of  influencing  the  performance  may  be  inves¬ 
tigated.  Close  tolerances  on  each  factor  must  be  held.  Extreme 
care  in  selection  of  primers  must  be  exercised.  Temperatures  should 
be  closely  controlled  and  reproducible. 

D.  The  firing  program  should  be  conducted  using  a  gran  barrel  sufficiently 
long  to  insure  that  "all  burnt"  occurs  well  before  shot  ejection. 

E.  Pressure  data  for  several  points  in  the  chamber  and  barrel  should  be 
obtained. 

F.  The  experimental  data  should  be  thoroughly  analyzed  to  determine 
what  possible  burning  and  Ignition  theories  may  apply  and  how  they 
may  be  applied  to  the  new  data. 
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APPENDIX  B 

AN  IDEALIZED  PRESSURE  -TIME  CURVE 


Consider  a  pres sure -time  curve  which,  from  a  given  shot  start  pressure 

rises  linearly  with  time  and  then  remains  constant  after  attainment  of  a 
high  pressure. 

The  ass umed pres  sure  variation  is 


=  1  +  bt 


b  =  a  constant 

P0  =  shot  start  pressure 


(B-l) 


Equation  of  motior  is  written 


AP  =  m'  ^ 
at 


(B -2) 


Integrate  (B-2)  using  (B-l) 


jrg  [(bt)  +  |  (bt)2_| 


Integrate  (B-3)  to  yield  travel 


(B-3) 


AP, 


x-x  = 

o 


(bt)' 


m'b2  L—  + 
The  energy  and  state  equations  are 


(bt)3 


(B-4) 


NR  (Tq-T)  =  m,v2 


(B-5) 


P  [  AX  .  . 


NRT 


f  (B-6) 
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Eliminate  T  betvaan  (B-5)  and  (B. 6)  using  (B-l),  (B-2),  (b-3),  and  (B-h). 
For  the  high  volumetric  loading  denaltlea  used  In  liquid  propellant,  the 


approximation  is  made  that  AX  = 

o 


N  - 


, «—  A2p  2 
(y-i)  A  po 

A 

m'b  J 

bt  + 

L  2  J 

2  A%S 

+  - %-  (l+bt) 

m'b^  1 

(B-7) 

F  (1  +  f ) 

The  quantity  y  can  be  neglected.  For  representative  propellant  at  50,000 
BP 

psi,  y  -  0.06  and  is,  of  course,  smaller  at  lower  pressure. 

Equation  of  burning  is  rewritten 


df  *  J 


(B-8) 


Differentiate  (B-7),  neglecting  ^  ,  and  equate  to  (b-8). 


A2po 

m'Fb 


-  .2.2 

y  (bt  +  — )  + 


c- 


b2t2  bV  ^  - 

O  +  —  r 


(l+bt) 


') 


(B-9) 


If  the  starting  pressure  is  not  unreasonably  high  and  if  the  rate  of  pres¬ 
sure  rise  is  high,  as  is  the  case,  bt  »  l  away  from  the  origin. 

AJter  (B-9)  accordingly. 


BS?  A2po 
2  =  m'Fb 


[7f+|)  (bt)  +  (-^ti)  (bt)2J 


(B-10) 
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From  (B-3) 


Integrate  (B-2)  using  (B-l6) . 


v  -  vl  -  -W  (*-*1)  (8-1T) 

Subscript  1  is  evaluated  when  the  constant  pressure  epoch  is  first  achieved. 

Eliminate  temperature  from  (B-5)  and.  (B-6) .  Differentiate  with  respect 
to  time,  using  (B-l6),  (B-17),  (B-l8),  and  (B-2)  yielding  (B-19) . 


Equate  (B-19)  to  (B-8) 


V  . 

dt 


(B-18) 


dN 

dt 


7  ap  v 

aP 

F  (H-yi) 


(B-19) 


BSf 
2  = 


y  av 


ap 

F  CU-ji) 


(B-20) 


Combine  (B-20)  with  (b-14)  and  differentiate 


O'2  ( 


(B-21) 
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PERMISSIBLE  FORM  FACTORS 
FOR  THE  CUBIC  FORM  FUNCTION 

Consider  the  cubic  form  function 

0=1-  (l+k2+k3)f  +  kgf2  +  k^f3 

From  (B-l4)  and  (C-l) 

~2  *  U+k2+k  )  -  2k2f  -  3k  f 2 
C  ' 

2 

C  *  a  positive  constant. 

The  criterion  for  validity  of  the  form  factors  is 

S>0  for  0<f<l 
(C-2)  at  f  m  1  and  f  =  0,  using  (C-3). 

1  +  k2  +  k^  >  0 


(C-l) 


(C-2) 


(C-3) 


(C-4) 
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Equations  (C-6)  and  (C-7)  are  necessary  but  not  sufficient  conditions 
for  validity  of  the  form  factors. 


The  geometrical  approach  can  be  used.  If  (C-6)  and  (C-7)  are  valid,  the 
surface  is  non-negative  at  the  end  points  of  the  interval. 

Equation  (C-2)  has  either  one  maximum  or  one  minimum  and  monotonically 
varies  from  the  zero  slope  position.  Therefore,  knowing  the  surface  at  the 

end  points  and  the  surface  and  location  of  the  zero  slope  position,  (C-3) 
can  be  checked. 

The  zero  slope  position  of  the  quadratic  S  curve  and  the  surface  there 
are  found  from  (C-2).  They  are  denoted  by  asterisks. 


f *  = 


(C-8) 


C2  '  1  +  k2  +  k3  +  f 


(C-9) 


*2  (A) 

c" 


(C-10) 


If  k3  >  0,  a  maximum  behavior  is  indicated  and  therefore,  the  validity 
of  (C-6) ,  (C-7)  and  (C-ll)  show  that  the  constants  are  valid. 


k^  >  0 


(C-ll) 
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If  k3  -  0,  the  cubic  degenerates  into  a  quadratic  which  is  of  no  in- 
teres t  here. 

If  k3  <  0,  a  minimum  behavior  is  exhibited.  From  (C-8),  since  k  <  0, 
k2  <  0  places  the  minimum  at  f  <  0,  the  proving  validity  of  the  constLts. 
-k2 

If  3k3  >  ^  the  lies  to  the  right  of  the  interval  and  the  con¬ 

stants  are  proven.  If  0  <  f*  <  1,  (C-9)  must  be  evaluated.  If  S*  >  0, 
the  form  factors  are  valid;  if  S*  <  0,  the  fonn  factors  are  invalid. 

— lng  UP-  lf  ValU6S  °f  k2  601  k3  are  examined  for  validity,  proceed 
as  follows; 

1.  Test  (C-6)  and  (C-?) 

a.  If  either  is  Invalid,  reject  the  form  factors. 

b.  If  both  are  valid,  proceed  to  step  2. 

2.  Test  (C-ll). 

a.  If  valid,  the  form  factors  are  proven. 

b*  If  k3  =  °'  the  foTa  function  is  quadratic  and  the  solution 
is  known. 

c.  If  Invalid,  proceed  to  step  3. 

3-  If  k2  <  0,  the  form  factors  are  proven 
k.  Evaluate  f*  from  (C-8) 

a.  If  f*  >  1,  the  form  factors  are  proven 

b.  If  f*  <  1,  proceed  to  step  5. 

5-  Compute  S*/c 2  from  (C-9) 

a.  If  S*  >  0,  the  form  factors  are  proven 

b.  If  S*  <  0,  the  form  factors  are  rejected 
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APPENDIX  d 


SOLUTION  OF  BALLISTIC  EQUATION 
USING  A  CUBIC  FORM  FUNCTION 

The  working  ballistic  equation  is  written 


-.a 

and 

(I2) 

(D-l) 

n 

1 

9 

M 

(D-2) 

$  =  1  -  (l+k2+k3)f  +  k2f2  +  k3f3 

(D-3) 

ii 

> 

x 

M 

(D-4) 

&t  -  —  =  a! 

\p  21 

(D-5) 

^3  =  form  factors 

The  above  ballietle  equations  are  subject  to  initial  conditions. 

1  '  0  t.tB  t.fa 

f-fo  A-  1  *.0 

(D-6) 

Combine  (D-2)  with  (D-3) 

0  =  0o  +  k5  (o'!)  +  kg  (csl)2  -  k3  (<M)3 

(D-7) 
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Equation  (D-l6)  is  a  well  known  form  end  is  integrated  from  the  initial 
condition.  The  variable  Z  is  eliminated  using  (D-13). 


X  “  tst  ~  a  0 


LA‘  J~ 


,1  -  4s  -  a  A  0 
L  f  °  ° 


i 

J  -  *2  (I)  Je  HI  (D-17) 


Je  (I)  -  exp  J~  ^  u  (d-18) 

Equation  (D-17),  with  <f>  expressed  in  terms  of  I  as  given  by  (D-7),  is 
the  desired  velocity-travel  curve.  The  pressure  is  found  by  combining 
(D-l)  and  (D-4)  to  yield  (D-19) 


r 


g  0*  -  | 

o* 


(D-19) 


The  time  may  he  found  by  graphically  or  numerically  integrating  (d-5). 

As  an  aid  to  computation,  (D-lB)  Is  integrated.  ra.  pertlal  fraetlon 
expansion  of  the  integrand  of  (D-l8)  may  be  found. 


*?*r-rT ij 


(D-20) 


ldiere  I  ,  j  and  I  are  the  3  roots  of  v  ..  , 

12  3  0  00‘'s  or  F1  (I),  which  can  be  found  by  well 

known  methods 


Al  = 


-  I, 


^  0s  Ix  -  3^3  o*  I12 


(D-21) 
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From  (D-l8)  and  (D-20) 


The  integral  in  (D-17)  can  be  evaluated  numeric®'* 

.j 

the  procedure  evaluating  (D-2U)  is  considered  too  te# 
numerically  integrated  for  Jc . 

*»• 

i. 


